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Abstract 

It is shown that Einstein- Yang-Mills-dilaton theory has a countable family of static globally 
regular solutions which are purely magnetic but uncharged. The discrete spectrum of masses 
of these solutions is bounded from above by the mass of extremal Gibbons-Maeda solution. As 
follows from linear stability analysis all solutions are unstable. 
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1 Introduction 

In a seminal paper [1] Bartnik and Mckinnon (BM) have discovered a countable family of glob- 
ally regular static spherically symmetric solutions of the Einstein- Yang- Mills (EYM) equations. 
A rigorous proof of existence of BM solutions was given by Smoller et al. [2]. Recently, Sudarsky 
and Wald have proposed a heuristic argument [3] , in the spirit of Morse theory, which explains 
existence and properties of BM solutions. This argument exploits the existence of topologi- 
cally nonequivalent multiple vacua in the SU(2)-YM theory (which is related to the fact that 
iT3(SU(2)) ~ Z) and is essentially insensitive to the concrete form of the coupling, which sug- 
gests that there should exist solutions similar to BM solutions in other theories involving the 
£?7(2)-YM field. Indeed, such solutions were found in YM-dilaton theory [4,5] and remarkable 
parallels between these solutions and BM solutions were observed [5]. 

The YM-dilaton theory and the Einstein- YM theory may be embedded in a single Einstein- 
YM-dilaton theory governed by the action 



where R is a scalar curvature, ^ is a dilaton, J 7 is a Yang-Mills curvature, G is Newton's 
constant and a is the dilaton coupling constant. This theory is characterized by a dimensionless 
parameter a = G/a 2 . When a = the action (1) reduces to the YM-dilaton theory. When 
a — > oo the action (1) becomes the Einstein- YM theory (plus trivial kinetic term for the scalar 
field). Finally, the case a = 1 corresponds to the low-energy action of heterotic string theory. 

The aim of this paper is to show that the theory defined by the action (1) has (for all values 
of a) static spherically symmetric globally regular solutions with the following properties: 

a) there exist a countable family of solutions X n (n £ N), 

b) the total mass M n increases with n and is bounded from above, 

c) the solution X n has exactly n unstable modes. 

The solutions X n depend continously on a and interpolate smoothly between the YM-dilaton 
solutions (a = 0) and BM solutions (a = oo). In the limit n — ► oo the solution X n {a) tends to 
the abelian extremal charged dilatonic black hole [6,7], whose mass therefore provides un upper 
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bound for the spectrum M n . These results were previously announced by the author in [5] and 
later the case a = 1 was analysed in [8]. 

The paper is organized as follows. In the next Section the field equations are derived and 
some scaling properties are discussed. In Section 3 the explicit abelian solutions are described. 
In Section 4 the numerical non-abelian solutions are presented and their qualitative properties 
are discussed. Section 5 is devoted to linear stability analysis. 

2 Field equations 

We are interested in static spherically symmetric configurations. It is convenient to parametrize 
the metric in the following way 



where A and N are functions of r. 

I assume that the electric part of the YM field vanishes (actually this is not a restriction 
because one can show, following the argument given in [9], that there are no globally regular 
static solutions with nonzero electric field). The purely magnetic static spherically symmetric 
SU(2) YM connection can be written, in the abelian gauge, as [10] 



where (i = 1,2,3) are Pauli matrices and w is a function of r. The corresponding YM 
curvature T = dA + eA A A is given by 



where prime denotes derivative with respect to r. 

For these anzatze and for (f> = <p(r) the action (1) gives the lagrangian (where S = 167T / Lett) 

L= / A(m'-U)dr , (5) 
Jo 

where the mass function m(r) is defined by N = 1 — 2Gm/r and 



ds 2 = -A 2 Ndt 2 + N^dr 2 + r 2 (cM 2 + sin 2 $d<p 2 ) , 



(2) 



eA = WTidfi + (cot $T3 + WT2) sin -Qd^p , 



(3) 



eT = w'ridr A dd + w'r2dr A siwddip — (1 — w^T^di) A siniMc^ 



(4) 



U = -r 2 N(t)' 2 + 4re" 2 ^ Nw' 2 + 




(6) 
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Hereafter, for convenience, I put the coupling constants e = a = 1 which means that I choose 
a/e as the unit of length and 1 /ea as the unit of energy. Then the system is characterized by 
one dimensionless parameter a = G/a 2 . Variation of L with respect of m, A, w, and (ft yields 
the field equations[| 



in 



aU , 



A 1 = aA(r(ft' 2 + -e 
r 



(ANe- 24> w')' + -^Ae- 2 *w(l 



w 



2 AN (/>')' + 2Ae~ 2 ^ 



2\2 



Nw 



,2 , (1-^) 



2r 2 



0. 



(7) 
(8) 
(9) 

(10) 



Note that L has a characteristic for general relativity "pure constraint" form, namely the 
integrand in Eq.(5) is the hamiltonian constraint, Eq.(7), multiplied by the lapse function A. 
Actually L is not differentiable because the variation of L gives an unwanted surface term at in- 
finity. To remedy this one has to add to L the Regge-Teitelboim correction term, — A(oo)m(oo). 
Then L' = L — A(oo)m(oo) has a well-defined functional derivative (since I deal with globally 
regular solutions I assume that m(0) = 0). It is convenient to define the energy functional 



E 



-L' 



{AU - A'm)dr , 



(11) 



which on shell is equal to the total mass m(oo) (assuming the boundary condition A{oo) = 1 
i.e. that t is proper time at spatial infinity). 

Eqs.(7)-(10) have a scaling symmetry: if w, (ft, N and A are solutions so are 



w\(r) = w(e x r) , 

<f> x (r) = (ft(e x r) + A , 
jV A ( r ) = N(e x r) m\(r) = e~~ x m(e x r) 
A x (r) = A(e x r). 

Under this transformation the energy functional scales as follows 



Ex 



(12a) 

(126) 
(12c) 
(12d) 

(13) 



The principle of symmetric criticality (see, R.Palais, Comm. Math. Phys. 69 (1979) 19) guaranties that the 
variation of S within the spherically symmetric ansatz gives the correct equations of motion. 
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The existence of this scaling symmetry does not exclude globally regular solutions because for 
the variation induced by the transformation (12) the energy is not extremized since 6(f)(oo) is 
nonzero. Hereafter, I will assume that all solutions satisfy 4>(oo) = 0, which can always be set 
by the transformation (12). This choice sets the scale of energy in the theory. 

3 Abelian solutions 

In the U{1) sector of the theory (i.e. in Einstein-Maxwell-dilaton theory) there are no static 
globally regular solutions but there are known explicit black hole solutions of Eqs.(7)-(10). The 
uncharged solution is Schwarzschild 

w = ±1 , 4> = , m = M = const , A = 1. (14) 

Charged black hole solutions were found by Gibbons and Maeda [6] (and later rediscovered in 
[7] ) . In the so called extremal limit the area of the horizon of these charged dilatonic black holes 
goes to zero and the resulting spacetime has a null naked singularity. This singular extremal 
solution will play an important role in our discussion of non-abelian solutions. It has a very 
simple form in isotropic coordinates 

ds 2 = -e'^dt 2 + e 2a<t> {dp 2 + p 2 d& 2 + p 2 sin 2 tidip 2 ) (15a) 

where 

and w = 0, so the YM curvature is 

T = -r 3 d-d A sin My , (15c) 

which corresponds to the Dirac magnetic monopole with unit magnetic charge (where the unit 
of charge is 1/e). The total mass is m(oo) = + a. For a = the solution (15) reduces 

to the dilatonic magnetic monopole discussed in [5]. For a = oo the solution (15) reduces to 
the extremal Reissner-Nordstrom black hole (before taking the limit a —> oo one has to make 
rescaling p — ► \fa,p). 
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4 Non-abelian solutions 

In order to construct solutions which are globally regular we have to impose the boundary 
conditions which ensure regularity at r = and asymptotic flatness. The asymptotic solutions 
of Eqs.(7)-(10) satisfying these requirements are 

± w = 1 - br 2 + 0(r 4 ) , (16a) 

<j> = c- 2e~ 2c b 2 r 2 + 0(r 4 ) , (166) 

N = 1 - 4ae" 2c 6V + 0(r 4 ) , (16c) 

A = d (1 + 4ae- 2c &V) + 0(r 4 ) (16d) 

near r = 0, and 

±™ = 1_^ + (1), (17a) 

<A = C-^+0(l), (176) 

Ar 1 2aM aL> 2 1. . . 

N = l - — + — + °i-z)i ( 17c ) 

^ = 1 -^ + 0(^) (17d) 

near r = oo. Here 6, c, d, i?, C, D, and M are arbitrary constants. All higher order terms in the 
above expansions are uniquely determined, through recurrence relations, by b, c and d in (16), 
and B,C,D and M in (17). The constant D is usually referred to as the dilaton charge while 
M = m(oo) is the total mass. 

Remark 1. Notice that, for the asymptotic behaviour (17a), the radial magnetic curvature, 
B r = ts(1 — w 2 ) / r 2 , falls-off as 1/r 3 , and therefore all globally regular solutions have zero YM 
magnetic charge. In accord with this, in the asymptotic expansion of 500 = ^4 2 iV the 1/r 2 term 
vanishes. 

Remark 2. One can easily show (see [5]) that for global solutions satisfying the above bound- 
ary conditions, a function w oscillates around zero between —1 and 1, while (j) is monotonically 
decreasing. 

Now, let us assume that near r = there exist a family of local solutions defined by the 
expansion (16). Note that this is a nontrivial statement because the point r = is a singular 
point of Eqs.(7)-(10), hence the formal power-series expansion (16) may have, in principle, a 
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zero radius of convergence. The initial parameters c and d are irrelevant since they can be 
chosen arbitrarily by the scaling (12) and by time rescaling, respectively. Thus, effectively we 
have a one-parameter family parametrized by b. For generic b the solution will not satisfy 
the asymptotic conditions (17) (in fact, the solution may even become singular at some finite 
distance). The standard numerical strategy, called the shooting method, is to find an initial 
value b for which the local solution extends to a global solution with the asymptotic behaviour 
(17). I have found that for generic orbits with b < &oo(a) a function w oscillates finite number 
of times in the region between w = — 1 and w = 1 and then goes to ±oo. For b > b OQ (a) all 
orbits become singular at a finite distance (in a sense that w' becomes infinite). 

The numerical results strongly indicate that for all values of a there exist a countable family 
of initial values b n (n £ N) determining globally regular solutions. Here the index n labels the 
number of nodes of the function w. In Table 1 are displayed the initial values and masses of the 
first five solutions for a = 1. The initial value c is determined by the condition that the dilaton 
vanishes at infinity i.e. C = 0. The functions w, 4> and N are graphed in Figs. 1-3. 

Table 1: Initial data (b,c) and masses of the first five globally regular solutions for a = 1. 



n 


be- c 


c 


M 


1 


0.1666666 


0.9311 


0.5773 


2 


0.2318001 


1.7925 


0.6849 


3 


0.246861 


2.6919 


0.7035 


4 


0.249483 


3.5974 


0.7065 


5 


0.249915 


4.5043 


0.7070 



In Figs. 4 and 5 are shown the functions w\ and N\ for three different values of a. For all 
a the solutions display three characteristic regions. The energy density is concentrated in the 
inner core region r < R\, where R\ is approximately the location of the first zero of w. In the 
second region, R\ < r < R2, where R2 is approximately the location of the last but one zero 
of w, the function w slowly oscillates around w = 0. Finally, in the asymptotic region r > R2, 
the function w goes monotonically to w = ±1 (hence the YM magnetic charge is gradually 
screened) and for r — > 00 the solution tends to the Schwarzschild solution (14). 
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Let us consider the limit n — > oo. In this limit R± — > and i?2 — ► oo, so the second region 
covers the whole space. As n grows the amplitude of oscillations of the function w decreases 
and goes to zero as n — ► oo (see Fig.l). Thus, for n — > oo the solution tends (nonuniformly) to 
the extremal Gibbons-Maeda solution (15). This is clear from Figs. 2 and 3, where the solution 
(15) for a = 1, expressed in Schwarzschild coordinates 



is included for comparison with non-abelian solutions. 

For all n and a the metric coefficient N has one minimum N m i n located approximately at 
R\. For fixed n, N m i n decreases with a (see Fig.5). Also for given a, N m i n decreases with n (see 
Fig.3) and N m i n — ► 1/(1 + a) 2 as n — > oo. 

The total mass M ra increases with n and for n — ► oo goes to l/\/l + a. Because our scale 
of energy breaks down for a — > oo (i.e. a = 0), it is more convenient to compare masses of 
solutions with different values of a using 1/eV G + a 2 as the unit of energy instead of 1/ea. 
This corresponds to rescaling M = \f\ + a M. Then for all a one has lirrin^ooMn = 1. It turns 
out that for given n the mass M n is increasing approximately linearly with (3 = a /(I + a). 
This is shown in Fig.6. In [5] it was proven in the case a = that the energy is equal to minus 
dilaton charge D. This property still holds for a > 0. To show this I first derive a simple virial 
identity. Consider a one-parameter family of field configurations defined by 




(18) 



and 




(19) 



w x (r) = w(Xr) , 



Mr) = <XAr) , 



N x (r) = N(Xr) 



A x {r) = A(Xr) . 



For this family the energy functional (9) is 



E x = A ii + A I 2 



(20) 



where 




(21) 




dr . 



(22) 
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Since the energy functonal is extremized at A = 1, it follows from (20) that on shell 



h = h- 



(23) 



Integrating Eq.(10) one gets D = —2/2, and therefore Eq.(23) yields 



to(oo) + D = . 



(24) 



Thus the total mass can be read off from the monopole term in the asymptotic expansion (17b) 
of the dilaton field. The identity (24) is also true for the limiting solution (15). 

5 Stability analysis 

In this Section I consider the time evolution of linear perturbations about the static solutions 
described above. I will assume that the time-dependent solutions remain spherically symmetric 
and the YM field stays within the ansatz (3). This is sufficient to demonstrate instability 
because unstable modes appear already in this class of perturbations. For radial perturbations 
the stability analysis is relatively simple, because the spherically symmetric gravitational field 
has no dynamical degrees of freedom and therefore the perturbations of metric coefficients are 
determined by the perturbations of matter fields [11]. To see this consider the Einstein equations 



where the functions v and A are defined by e" = AN and e A = N, and dot denotes the time 
derivative. Now, take the perturbed fields v{r) + 5v(r,t), A(r) + SX(r, t), w(r) + 5w(r,t), and 
<j>{r) + Scf)(r,t), where {y{r), A(r), w(r), <t>{r)) is a static non-abelian solution, and insert them 
into Eqs.(25) and (26). Linearizing Eq.(25) one obtains 



A = —2ar(4>4> H — ^ww ) 



(25) 



and 



v < - \< = ar U' 2 + %w' 2 + e- 2 »(4> 2 + ^w 2 ) 



(26) 



<5A 




(27) 



which can be integrated to give 



<5A 



2ar((f)'5(j) H — ^w'5w) . 



(28) 
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Linearization of Eq.(26) yields 



Si/ - 5X' = 2ar(<p'5<f>' + ^ w'Sw') . 



{r<f>" + <j>')5<i> + \w" - ^w')5w 



and therefore using Eq.(28) one gets 

5u' = -2a 

The spherically symmetric evolution equations for the matter fields are 

-{e- u - 2<t> w) + (e"- 2 *u/)' + ^e^-^wil - w 2 ) = , 

and 

-r 2 (e^>)-+ (r 2 e u 4>y + W^' 2 * 



2^2 



X ,2 {l-W 2 ) 



2r 2 



= . 



(29) 



(30) 



(31) 



(32) 



Multiplying Eqs.(31) and (32) by e linearizing and assuming harmonic time dependence for 
the perturbations, 5w(r, t) = e' iat ^(r) and 5<p(r, t) = e tat ip(r), one obtains an eigenvalue problem 



+ (2(1)' - + 2v/fJ/ + 2aw' 
2a 



(r<f/' + <t>'U + -(w" --v/)Z 
r r 

2 



1 



^e" A (l-3«; 2 )e 



+— e- A w(l - w 2 )(<f>'4> + -W£) = a 2 e- lv i , 



(33) 



_( r 2^)'_4 e -20 



-rV-0' + 2arV 



(r0" + 0')V' + ^K-^Oe 



2 2 -2i> / 

are tp . 



(34) 



If the perturbations satisfy the boundary conditions £(0) = 0, -0(0) = const, £(oo) = 0, and 
ip(oo) = 0, then the above system is self-adjoint, hence eigenvalues a 2 are real. Instability 
manifests itself in the presence of at least one negative eigenvalue. 

I have used the generalized rule of nodes (Jacobi criterion) to find that for all values of a 
the nth solution has exactly n negative eigenvalues (I have checked this up to n = 4). This is 
consistent with the fact that the limiting solution (15) has infinitely many unstable modes which 
can be seen easily by inserting (15) into Eq.(33) and considering perturbations with ip = 0. The 
result that the nth solution has n unstable modes is essential for the Sudarsky-Wald argument 
[3]. However, by considering a restricted class of perturbations some directions of instability 
might have been suppressed. 

Acknowledgement I am grateful to Peter Aichelburg for stimulating discussions. 
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Figure captions 

Fig.l The YM potential w n for n = 1,3, 5. 

Fig.2 The dilaton <f) n for n = 1, 3, 5. 

Fig.3 The metric coefficient for n = 1, 3, 5. 

Fig.4 w\ for a = 0.01 (dashed line), a = 1 (solid line) and a = 100 (dotted line). 
Fig.5 N\ for a = 0.01 (dashed line), a = 1 (solid line) and a = 100 (dotted line). 
Fig.6 The total mass M n (in units 1/ey/G + a 2 ) as a function of (3 = a/(l + a). 



